function  interp_function  = Interp_function_simulation(estim_value,State_vars)

log_Skills_tp1 = log(State_vars(:,1));
H_prime        = log(State_vars(:,2)); 

interp_function_hat(:,1)= estim_value(1).*ones(size(log_Skills_tp1));
interp_function_hat(:,2)= estim_value(2).*log_Skills_tp1;
interp_function_hat(:,3)= estim_value(3).*H_prime;
interp_function_hat(:,4)= estim_value(4).*log_Skills_tp1.*H_prime;
                             
interp_function = sum(interp_function_hat,2);      